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Abstract 

This  research  effort  investigated  characterizing  the  point  spread  function  (PSF) 
behavior  of  radially  displaced  point  scatterers  using  circular  synthetic  aperture  radar 
(CSAR).  For  persistent  staring  radar  applications  utilizing  CSAR  for  target  recog¬ 
nition,  it  is  important  to  know  how  a  target’s  PSF  behaves  as  a  function  of  various 
radar  functional  parameters  and  different  target  positions.  Thus  far,  research  has 
been  conducted  to  understand  PSF  of  a  scatterer  located  at  the  imaging  scene  center. 
An  analytic  closed-form  solution  has  been  derived  assuming  the  scatterer  is  located 
at  the  origin  of  the  CSAR  imaging  geometry.  However,  it  is  difficult  to  derive  an 
analytic  PSF  solution  for  a  scatterer  that  is  radially  displaced  from  the  imaging  scene 
center.  Using  the  backprojection  image  formation  algorithm,  PSF  responses  are  gen¬ 
erated  at  various  point  target  locations.  Consistent  with  previous  studies,  the  three 
dimensional  PSF  for  a  point  target  located  at  the  image  center  is  cone  shaped  and 
serves  as  the  basis  for  comparing  and  characterizing  the  PSFs  of  radially  displaced 
scatterers.  Simulated  results  show  the  impulse  response  of  a  radially  displaced  point 
scatterer  is  asymmetric  and  tends  to  exhibit  increased  cllipticity  as  it  moves  further 
from  the  scene  center. 
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Point  Spread  Function  Characterization 
of  a  Radially  Displaced  Scatterer 
Using  Circular  Synthetic  Aperture  Radar 

I.  Introduction 

Synthetic  Aperture  Radar  (SAR)  is  a  revolutionary  technology  invented  by  Carl  Wiley 
in  1950.  SAR  applications  benefit  both  the  civilian  and  military  sectors.  Examples 
of  SAR  applications  include  environmental  monitoring  [16],  terrain  mapping,  and 
surveillance  operations.  The  most  important  feature  of  SAR  is  its  capability  to  form 
images  in  day  or  night  under  all  weather  conditions.  For  staring  radar  applications, 
the  use  of  SAR  technology  is  paramount. 

Circular  Synthetic  Aperture  Radar  (CSAR)  is  a  method  for  collecting  raw  SAR 
data.  In  CSAR  applications,  the  radar  moves  in  a  circular  flight  path  while  keeping 
its  mainbeam  focused  at  the  center  of  the  ground  patch  to  be  imaged.  The  CSAR 
process  is  very  similar  to  inverse  SAR  (ISAR)  where  the  target  rotates  and  the  radar 
is  stationary.  Other  methods  for  collecting  raw  SAR  data  include  stripmap  SAR, 
spotlight  SAR  and  interferometric  SAR  (IFSAR).  In  stripmap  SAR  mode,  the  radar 
moves  along  a  straight  flight  path  and  illuminates  a  rectangular  strip  of  the  ground 
during  data  collection.  In  spotlight  SAR  mode,  the  radar  moves  along  a  straight  flight 
path  while  steering  its  mainbeam  to  constantly  illuminate  a  specific  ground  patch. 
Lastly,  IFSAR  mode  operates  by  using  data  collected  from  two  separate  receiver 
channels  that  are  vertically  displaced  on  the  platform.  The  height  differential  is  then 
used  to  map  the  imaging  scene  in  3D  using  interferometry. 

Each  of  the  SAR  modes  described  present  different  operational  issues  and  lim¬ 
itations  when  forming  the  desired  image.  For  example,  a  stripmap  image  is  formed 
using  data  collected  over  limited  azimuth  look  angles.  Thus,  stripmap  images  may  not 
be  of  sufficient  quality  to  permit  target  identification  with  a  high  degree  of  confidence. 
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For  an  automatic  target  recognition  (ATR)  system  requiring  target  identification  with 
a  higher  degree  of  confidence,  CSAR  processing  represents  a  viable  alternative  given 
it  can  produce  images  with  resolution  less  than  a  wavelength,  depending  on  scatterer 
persistence  throughout  the  collection  interval.  Many  studies  have  been  conducted 
that  address  CSAR  image  reconstruction  [1,3,7,18],  with  [19]  describing  the  CSAR 
image  reconstruction  process  in  great  detail. 

1.1  Research  Problem 

Wide-angle  imaging  provides  significantly  better  resolution  than  other  imaging 
modalities.  Higher  resolution  has  proven  to  provide  better  target  recognition  per¬ 
formance  for  ATR  systems.  For  that  reason,  CSAR  has  been  used  for  staring  radar 
applications  in  which  large  apertures  could  be  utilized  for  achieving  better  cross-range 
resolution.  In  previous  studies  considering  wide-angle  imaging  and  point  spread  func¬ 
tions  (PSF),  the  studies  assumed  that  the  target  was  located  at  the  imaging  scene 
center.  Based  on  this  assumption  and  the  assumption  that  the  entire  circular  path 
is  flown,  the  behavior  of  a  point  target  PSF  is  well  understood.  Researchers  have 
derived  analytic  closed-form  PSF  solutions  for  a  point  target  located  at  the  center  of 
the  imaging  geometry  [13,19]. 

One  important  issue  in  CSAR  imaging  involves  understanding  target  PSF  be¬ 
havior  as  the  target  becomes  radially  displaced  from  the  image  scene  center.  Complete 
PSF  characterization  of  a  specific  target  (i.e. ,  one  comprised  of  multiple,  uniquely  lo¬ 
cated  point  scattering  centers)  could  lead  to  enhanced  focusing  and  geolocation  of 
this  particular  target  as  it  moves  through  the  CSAR  imaging  geometry. 

It  is  well  understood  that  not  all  target  scattering  is  point-like  and,  therefore, 
additional  challenges  concerning  complex  target  characteristics  (number  of  dominant 
scatterers,  scattering  center  location,  orientation,  etc.)  should  be  investigated;  how¬ 
ever,  these  challenges  are  beyond  the  scope  of  this  research.  Thus,  this  research 
focused  on  characterizing  the  PSF  behavior  of  a  single  point  scatterer  that  is  first 
located  at  the  image  center  and  then  radially  displaced  within  the  scene. 
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1.2  Research  Motivation 


In  today’s  world,  the  battlefield  is  becoming  more  and  more  complex.  This  com¬ 
plexity  grows  as  battles  move  from  open  field,  to  forest,  to  urban  areas.  Designating 
targets  as  belonging  to  enemy  forces  is  becoming  more  and  more  difficult  as  well, 
given  that  military  vehicles,  civilian  vehicles,  bicycles,  etc.,  are  now  commonly  used 
to  conduct  military  operations.  Thus,  identification  of  complex  targets  in  complex 
battle  areas  is  very  important  and  extremely  challenging.  Persistent  staring  radar 
utilizing  CSAR  represents  one  promising  technology  to  address  the  aforementioned 
issues. 

Synthetic  Aperture  Radar  (SAR)  is  commonly  used  for  imaging  and  automatic 
target  recognition  (ATR)  applications.  CSAR  collections  with  large  apertures  (those 
spanning  90°  of  azimuth  or  more)  have  two  properties  that  potentially  make  them 
very  useful.  First,  they  can  provide  resolution  in  the  X-Y  planes  that  are  less  than  a 
wavelength  and  are  essentially  independent  of  the  bandwidth.  Second,  the  resultant 
PSFs  can  provide  resolution  in  three  dimensions  using  a  single  phase  center  and  a 
single  pass.  It  is  envisioned  that  a  better  understanding  of  PSF  characteristics  for 
a  radially  displaced  target  could  enable  an  ATR  system  to  fully  exploit  these  two 
properties. 

One  significant  challenge  for  CSAR  processing  is  the  imaging  of  large  scene  sizes. 
For  large  scenes,  the  PSF  response  can  change  significantly  as  the  point  of  interest 
is  radially  displaced  from  scene  center.  Thus,  it  is  important  to  understand  PSF 
behavior  for  all  locations  in  the  scene,  and  in  particular,  when  imaging  large  scenes. 

The  fundamental  goal  of  this  thesis  research  is  to  characterize  the  PSF  response 
of  a  radially  displaced  point  target.  The  research  utilized  simulation  to  perform  this 
characterization.  The  ability  to  understand  PSF  behavior  for  a  target  located  in  any 
position  of  a  large  scene  is  fundamental  to  understanding  the  signature  of  complex 
targets.  This  understanding  will  lead  to  a  better  ATR  algorithm  development  for 
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both  target  recognition  and  detection,  e.g.,  ground  moving  target  indication(GMTI) 
development. 

1.3  Thesis  Organization 

The  rest  of  this  document  is  divided  into  five  chapters,  with  each  providing 
specific  details  on  a  topic  contributing  to  the  research.  The  topics  addressed  include: 
SAR  background,  CSAR  processing,  data  generation,  results  and  analysis. 

Chapter  2  describes  SAR  image  processing  concepts  to  include  a  definition  of 
range  and  cross-range  resolution.  Range  and  cross-range  imaging  is  outlined.  SAR 
data  collection  modalities  such  as  stripmap  SAR,  spotlight  SAR,  interferometric  SAR 
(IFSAR)  have  been  discussed  in  great  details. 

Chapter  3  provides  a  comprehensive  discussion  on  CSAR  processing.  Details 
on  an  important  SAR  image  formation  algorithm,  known  as  convolutional  backpro- 
jection,  are  presented  in  this  chapter. 

Chapter  4  describes  the  data  generation  methodology  used  to  characterize  the 
PSF  of  radially  displaced  point  scatterers.  CSAR  data  collection  geometry  is  intro¬ 
duced  and  key  CSAR  parameters  for  synthetic  data  generation  are  presented,  includ¬ 
ing  depression  angle,  radar  center  frequency,  bandwidth,  range  from  target  to  the 
aircraft.  Radar  waveform  selection  issues  are  also  discussed. 

Chapter  5  provides  simulation  results  and  analysis  for  several  scenarios.  A 
baseline  PSF  response  is  first  considered  where  the  point  target  is  located  at  the 
scene  center.  Both  2D  and  3D  PSFs  are  presented  for  the  baseline  case.  Results  for 
scenarios  involving  radially  displaced  point  targets  are  then  introduced.  Again,  both 
2D  and  3D  PSF  responses  are  presented  and  discussed. 

Chapter  6  provides  research  conclusions  and  presents  recommendations  for  fu¬ 
ture  research. 
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II.  Background  on  Synthetic  Aperture  Radar 


This  section  describes  synthetic  aperture  radar  (SAR)  signal  and  image  processing 
concepts.  The  SAR  signal  processing  concept  is  based  upon  single  (e.g.  range 
profile)  or  multi  dimensional  radar  signal  processing.  Many  books  have  been  written 
that  explain  radar  signal  processing  concepts  in  great  details.  The  radar  signal  pro¬ 
cessing  concepts  presented  in  this  chapter  are  extracted  from  books  written  by  Merrill 
Skolnik  [17],  Mark  Richards  [14]  and  Bassern  Mahafza  [10].  Several  well  written  SAR 
monographs  have  been  authored  by  Mehrdad  Soumekh  [19],  Charles  Jakowatz  [9]  and 
Walter  Carrara  [4].  Soumekh  [19]  provides  a  very  comprehensive  signal  level  under¬ 
standing  of  SAR.  Hence,  the  notation,  symbols,  mathematical  equations  and  pictures 
used  in  this  section  to  illustrate  SAR  signal  and  image  processing  concepts  follow 
from  Soumekh’s  book  [19].  Since  this  thesis  research  focused  on  CSAR,  a  limited 
discussion  on  basic  radar  concept  has  been  provided. 

2.1  Basic  Radar  Concept 

RADAR  stands  for  RAdio  Detection  And  Ranging.  A  radar  system  transmits 
a  sequence  of  pulses  which  are  incident  upon  targets  of  interest  and  reflect  back.  A 
receiver  records  the  echoed  signal.  Time  delay  between  the  transmit  and  receive  signal 
is  used  to  calculate  range  information.  Radar  range  is  calculated  using  the  following 
equation  [10, 14, 17] 


R=~f’  (2-1) 

where  c  is  the  speed-of- light  and  Tr  is  round  trip  travel  time  of  a  radar  signal. 

Figure  2.1  shows  the  radar  signal  transmission  and  reception  concept  [21].  After 
signal  transmission,  the  radar  system  waits  for  sufficient  time  to  record  the  echo 
signal.  The  waiting  period  (time)  between  two  consecutive  transmitted  pulses  is  called 
pulse  repetition  interval  (PRI).  The  reciprocal  of  PRI  is  called  the  pulse  repetition 
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frequency  (PRF).  Multiple  PRIs  are  necessary  to  detect  two  targets  unambiguously. 
The  following  equation  defines  maximum  unambiguous  range  [17] 


Figure  2.1:  Radar  signal  transmission  and  reception  concept  [21]. 


A  series  of  echoed  radar  pulses  form  a  coherent  processing  interval  (CPI)  or 
dwell.  Each  pulse  within  a  CPI  is  sampled  adequately  to  retrieve  meaningful  target 
information.  Each  sampled  point  of  a  pulse  is  known  as  range  bin,  range  gate,  or 
resolution  cell  [14].  Since  a  radar  system  samples  a  pulse  as  it  arrives,  range  bin  is 
also  known  as  fast  time  sample  and  range  axis  is  also  known  as  fast  time  axis.  After 
sampling  the  first  pulse,  the  next  pulse  is  sampled  and  the  process  continues  until  the 
last  pulse  of  the  CPI  is  received.  Hence,  the  pulse  axis  is  known  as  slow  time  axis. 
In  SAR,  slow  time  is  also  known  as  synthetic  aperture  axis.  Figure  2.2  shows  a  2D 
radar  data  matrix. 

Multiple  receive  channels  provide  multiple  CPI  processing  and  constitute  a  data 
cube.  Figure  2.3  shows  a  3D  radar  data  cube.  The  3D  radar  data  cube  concept  is 
illustrated  by  Mark  Richards  [14]  and  provides  an  easy  way  of  understanding  and 
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Slow  time  or  Pulse  number 


Figure  2.2:  2D  Radar  Data  Matrix.  The  horizontal  axis  is  the  slow  time  or  pulse 
number  axis.  It  contains  M  pulses  for  a  CPI  and  represents  duration  of  CPI.  The 
vertical  axis  is  the  fast  time  or  range  axis.  It  contains  L  range  bins  or  fast  time 
samples  and  represents  duration  of  CPI  [14], 
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remembering  several  important  radar  signal  processing  concepts.  For  example,  ID 
filtering  or  Fourier  transform  along  pulse/slow  time  axis  is  called  Doppler  processing. 
Doppler  processing  provides  information  about  moving  targets.  Figure  2.4  shows 
Doppler  processing  concept. 

Pulse  compression  is  the  ID  convolution  along  the  range  axis  or  fast  time  axis. 
To  detect  targets  at  a  long  range,  very  long  pulses  are  needed,  ffowever,  a  long 
pulse  provides  poor  range  resolution.  To  mitigate  this  problem,  a  frequency  or  phase 
modulated  signal  is  transmitted  by  a  radar.  This  is  known  as  pulse  compression 
process.  Pulse  compression  enables  obtaining  the  high  range  resolution  of  a  short 
pulse  while  using  a  long  pulse. 

One  dimensional  (ID)  weighting  across  antenna  receive  channels  is  called  beam¬ 
forming.  Two  Dimensional  (2D)  joint  adapting  weighting  across  antenna  receive  chan¬ 
nels  and  pulse  numbers  is  called  space-time  adaptive  processing  (STAP).  Finally,  2D 
matched  filtering  in  slow  time  and  fast  time  is  called  synthetic  aperture  radar  (SAR) 
imaging. 

2.2  SAR  Imaging 

The  fundamental  SAR  imaging  geometry  is  shown  in  Figure  2.5  [15].  In  this 
figure,  the  aircraft  carrying  the  radar  flies  a  straight  flight  path  with  velocity  v  at 
a  certain  altitude.  Slant  range  is  the  distance  from  the  radar  to  the  imaging  scene 
center.  Azimuth  angle  0  is  defined  to  be  the  angle  around  the  scene  center  which 
encompasses  0  to  360°.  Azimuth  angle  is  also  known  as  aspect  angle.  Elevation  angle 
6  is  defined  to  be  the  angle  measured  from  the  ground  plane  toward  the  positive  z  - 
axis  which  encompasses  0  to  90°.  The  angle  between  the  slant  range  and  the  ground 
plane  is  a  typical  measure  of  the  elevation  angle.  Elevation  angle  is  equivalent  to  the 
depression  angle.  Ground  squint  angle  is  the  angle  between  the  aircraft’s  flight  path 
to  the  projection  of  the  slant  range  on  the  ground. 


Figure  2.3:  Radar  Data  Cube.  The  horizontal  axis  is  the  pulse  number  or  slow 

time  axis.  It  contains  M  pulses  for  a  CPI.  The  axis  into  the  page  is  the  fast  time  or 
range  axis.  It  contains  L  fast  time  samples.  The  vertical  axis  is  the  number  of  receive 
channel.  It  contains  N  receive  channels. 


L 


Apply  FFT  to  each  range  bin 
across  slow  time 

M 

1 

1  - ► 

Slow  time  or  Pulse  number  Frequency  bin  or  Doppler  bin 


Figure  2.4:  Range  Doppler  Processing  concept.  By  applying  Fourier  transform 

along  pulse/slow  time  axis  for  each  range  bin,  a  range  Doppler  map  is  created. 
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Figure  2.5:  SAR  imaging  geometry  and  terminology.  Azimuth  angle  cj)  is  the  angle 
around  the  scene  center.  Elevation  angle  9  is  defined  to  be  the  angle  measured  from 
the  ground  plane  to  the  z  -axis  [15]. 

2.2.1  Range  Imaging.  System  geometry  for  range  imaging  is  shown  in 
Figure  2.6  [19].  As  shown,  targets  are  located  at  X\,  X2  and  Xn.  Each  target  has  a 
reflectivity  of  a\,  a2  and  an.  Total  target  area  coverage  is  2X0.  For  range  imaging, 
all  targets  are  assumed  to  be  located  at  the  same  cross-range  location. 


Figure  2.6:  SAR  range  imaging  geometry.  Targets  reside  within  the  finite  spatial 

region  [Xc-Xq,  Xc+Xq]  or  range  domain.  This  is  also  known  as  the  range  swath.  Xc 
is  the  scene  center.  Total  target  spatial  support  is  2X0  [19]. 


Consider  p(t)  as  the  radar  transmit  signal.  Then  the  received  signal  s(t)  is 
characterized  by  the  following  equation: 
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s(t)  =  <JnP  ft  -  ^p) 


(2.3) 


where  xn  and  an  are  the  range  and  reflectivity  of  the  nth  target. 

Consider  S(u)  and  P(u)  are  the  Fourier  transform  of  s(t)  andp(f),  respectively. 
Therefore,  target  response  function  f(x)  can  be  reconstructed  using  the  following 
equation: 


/W  = 


P(u). 


(2.4) 


where  is  the  inverse  Fourier  transform  operator. 

Target  reconstruction  using  (2.4)  is  known  as  deconvolution.  This  method  is 
considered  numerically  unstable  when  P(uj)  is  zero.  For  that  reason,  matched  filtering 
is  used  for  target  reconstruction.  Also,  matched  filtering  provides  the  highest  signal 
to  noise  ratio  (SNR)  when  the  measured  signal  is  corrupted. 

Using  matched  filtering,  target  response  function  f(x)  is  given  by  [19] 


f(x )  =  ^  crnpsfn(x  -  xn) 

n 

where  psf(x)  is  defined  by  F]r'  [|P(o;)|2]. 


(2.5) 


2. 2. 1.1  Range  Resolution.  The  SAR  range  resolution,  Ax  is  given  by 
following  equations: 


Ax  = 


cAt 

2 


c 

2  u> 


(2.6) 


where  At  is  the  sampling  interval  or  PR1  and  to  is  the  bandwidth  in  Hertz. 

Figure  2.7  illustrates  the  range  image  reconstruction  via  matched  filtering.  The 
software  used  for  this  simulation  is  provided  by  Soumekh  [19].  Scene  center  Xc  is 
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located  at  the  range  2000  meters.  Targets  are  located  at  -25,  0,  17.5  and  35  meters 
from  the  scene  center  Xc  (i.e.  at  the  range  of  1975,  2000,  2017.5,  and  2035  meters). 
Reflectivity  of  these  targets  are  0.8,  1.0,  1.0  and  0.8,  respectively.  The  radar  center 
frequency  is  1.0  GHz  and  the  bandwidth  is  100  MHz.  As  shown  in  Figure  2.7,  target 
responses  represent  reflectivity  used  in  the  simulation. 


Figure  2.7:  Range  image  reconstruction  via  matched  filtering.  Scene  center  Xc  is 
located  at  range  2000  meters.  Targets  were  located  at  -25,  0,  17.5  and  35  meters  from 
the  scene  center. 

2.2.2  Cross-range  Imaging.  The  system  geometry  for  cross-range  imaging 
is  shown  in  Figure  2.8  [19].  In  cross-range  imaging,  the  radar  stays  at  a  fixed  range, 
x.  The  aircraft  moves  in  cross-range  axis,  y. 

As  depicted  in  Figure  2.8,  targets  are  located  at  Yl,  Y2  and  Yn.  Each  target 
has  a  reflectivity  of  ay,  a2,  and  an.  Total  target  area  coverage  is  2 Y0.  Synthetic 
aperture  axis  u  extends  from  (-L,  L).  Consider  the  case  when  the  radar  is  located  at 
(0,  u )  which  is  equivalent  to  ( x ,  y )  in  Cartesian  coordinate.  The  radar  distance,  d 
from  a  target  at  ( x ,  yn )  is  given  by  the  following  equation: 
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Figure  2.8:  Cross-range  imaging  geometry.  Targets  reside  within  the  finite  target 

area  [ — Vq,  Vq]  in  cross-range.  Total  target  cross-range  support  is  2 Y0  [19]. 


d  =  \/x2  +  (yn  -  u)2 


(2.7) 


Now  assume  the  radar  transmits  the  following  complex  signal: 


p(f)  =  exp  (jut) 


(2.8) 


The  measured  echo  signal  from  (2.8)  is  a  function  of  u  which  could  be  expressed  by 
the  following  equation: 


KC  u)  =  °nV 


t  — 


2^/x2  +  tjh  -  u)2 


(2.9) 


For  radar  frequency  u,  wave  number  k  is  defined  by  the  following  equation: 


c 


(2.10) 
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It  follows  from  (2.8)  and  (2.10)  that  (2.9)  can  be  rewritten  as  follows: 


s(t,  u) 
s[t,  u ) 


exp  {jut)  ^2  an  exp 

n 

exp  (  jut)  ^2  <rn  exp 

n 


0 .  Vx2  +  (yn~u)2 

-2  JUJ- - 

c 

-j]  2 k  y/x2  +  (yn  -  u)2 


After  baseband  conversion  using: 


(2.11) 

(2.12) 


s(u;,  u)  —  s(t,  u)  exp(—  jut)  (2-13) 

The  signal  in  (2.12)  can  be  rewritten  as: 

s{u,u)  =  ^  on  exp[— j'2fc\/ x2  +  (yn  -  u )2]  (2.14) 

n 

Now,  let  sn(u,u )  =  exp [—j'lk^Jx1  +  (yn  —  u)2].  Then  the  baseband  signal  in  (2.14) 
is  represented  by: 


s{u,u)  =  \  onsn{u,u) 


(2.15) 


The  Fourier  transform  of  the  phase-modulated  baseband  signal  in  (2.14)  is  given 
by  the  following  equation  [19]. 


Sn(uj,  ku)  =  on  exp  (-j  yj  Ak 2  -  klxn  -  jkuyn )  (2.16) 

where  ku  is  the  spatial  frequency  of  u  and  ku  E  [—2k,  2k].  ku  is  also  known  as  the 
slow-time  or  Doppler  domain  [19]. 

Let  s(ui,  ku)  be  the  sum  of  the  spatial  Fourier  transform  of  the  individual  signals 
sn(u>,ku).  In  cross-range  imaging,  all  the  targets  are  located  at  same  range  value. 
Therefore,  xn  =  Xc.  Hence  (2.16)  can  be  rewritten  as  follows: 
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(2.17) 


s(u,  ku)  =  unexp(-j^4k^k^Xc  -  jkuyn) 

It  can  be  shown  that  the  spatial  Fourier  transform  of  the  reference  signal  with 
respect  to  the  synthetic  aperture  u  is  the  following  [19]: 

s0(w,  ku)  =  exp(-jy/4^2  -  k\Xc )  (2.18) 

where  ku  e  [—2k,  2k]  and  zero  otherwise. 

Using  matched  filtering,  the  target  function  f(y)  in  cross-range  can  be  con¬ 
structed  per  the  followings: 


F(kt) 

=  s(u,ku)sl(u;,ku) 

(2.19) 

F(ku) 

=  ^2anexp(-jkuyn) 

(2.20) 

f(y) 

TL 

=  r^inK)] 

(2.21) 

where  F^k  1  ^  is  the  inverse  Fourier  transform  operator  and  Sq  is  the  complex  conjugate 

of  s0- 


f(y)  =  ^2  arlpsfn(y  -  yn)  (2.22) 

n 

Figure  2.9  shows  cross-range  image  reconstruction  via  matched  filtering.  Again, 
software  used  for  this  simulation  is  provided  by  Soumekh  [19].  Scene  center  Yc  is 
located  at  y  =  0  meters.  Targets  are  located  at  -50,  -25,  25  and  25  meters  from  the 
scene  center.  Reflectivity  of  these  targets  are  0.8,  1.0,  1.0  and  0.8,  respectively.  The 
radar  center  frequency  is  1.0  GHz  and  the  bandwidth  is  200  MHz.  As  shown  in  Figure 
2.9,  target  responses  represent  reflectivity  used  in  the  simulation. 
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Figure  2.9:  Cross-range  image  reconstruction  via  matched  filtering.  Scene  center 

Yc  is  located  at  y  —  0  meters.  Targets  are  located  at  -50,  -25,  25  and  25  meters  from 
the  scene  center. 

2.2.2. 1  Cross-range  Resolution.  The  SAR  cross-range  resolution  Ay 
is  given  by  following  equation  [19]: 

=  j  (2.23) 

where  A  =  A-  is  the  wavelength  of  radar  center  frequency. 

2.3  Stripmap  SAR 

Stripmap  SAR  was  the  first  SAR  imaging  method  invented  in  1950s  and  is  also 
known  as  side-looking  SAR.  In  stripmap  SAR  mode,  the  radar  moves  along  a  straight 
flight  path  and  illuminates  a  strip  of  the  ground  during  data  collection. 

Fundamental  Stripmap  imaging  geometry  is  shown  in  Figure  2.10  [19].  Radar 
carrying  the  aircraft  flies  along  the  cross-range  maintaining  a  fixed  range.  Data  is 
collected  and  image  is  formed.  Some  image  formation  algorithms  (e.g.  Fresnel  ap¬ 
proximation)  are  not  suitable  for  stripmap  SAR  processing  [19]. 
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For  some  applications  it  is  important  to  scan  some  areas  of  interest  quickly  to 
detect  target  presence.  Based  on  the  finding  a  target  is  present,  the  radar  operator 
may  employ  a  high  resolution  SAR  (e.g.  Spotlight  Mode  SAR)  system  for  obtaining 
more  detailed  target  information.  Stripmap  SAR  is  used  for  these  kind  of  applications. 


Figure  2.10:  Stripmap  data  collection  geometry.  Radar  carrying  the  aircraft  flies 

along  the  cross-range  maintaining  a  fixed  range.  Data  is  collected  and  image  is 
formed  [19]. 

2.f  Spotlight  SAR 

In  spotlight  SAR  mode,  the  radar  moves  along  a  straight  flight  path  while 
steering  its  mainbeam  to  constantly  illuminate  a  specific  ground  patch.  SAR  spotlight 
mode  imaging  geometry  is  shown  in  Figure  2.11  [19]. 

The  radar  carrying  aircraft  flies  in  the  cross-range  direction  (y-axis)  while  main¬ 
tains  the  fixed  range (x-axis).  Mechanical  or  electronic  beam  steering  is  used  to  focus 
radiation  pattern  on  the  same  target  area  [19].  Spotlight  SAR  data  collection  enables 
high  resolution  imaging  capabilities. 

2.5  Interferometric  SAR  (IFSAR) 

Interferometric  SAR  (IFSAR)  mode  operates  by  using  data  collected  from  two 
separate  receiver  channels  that  are  vertically  displaced.  The  height  differential  is  then 
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Figure  2.11:  Spotlight  SAR  data  collection  geometry.  The  radar  carrying  aircraft 
moves  along  a  straight  flight  path  while  steering  its  mainbeam  to  constantly  illuminate 
a  specific  ground  patch  [19]. 

used  to  map  the  imaging  scene  in  3D  using  interferometry.  Two  different  methods 
of  IFSAR  data  collection  are  found  in  literature.  Jakowatz  [9]  describes  these  two 
methods  as  single-pass  and  two-pass  interferometry. 

Single-pass  IFSAR  data  collection  geometry  is  shown  in  Figure  2.12.  In  this 
method,  two  receivers  reside  on  the  same  aircraft  with  a  small  vertical  separation. 
Two  images  can  be  formed  at  the  same  time  and  a  3D  terrain  map  can  be  generated. 
Single-pass  IFSAR  provides  better  registration  accuracy  for  interferometric  processing 
and  is  considered  to  be  very  useful  for  moving  target  indication  (MTI)  applications  [9] . 

Two-pass  IFSAR  data  collection  geometry  is  shown  in  Figure  2.13.  In  this  case, 
the  radar  first  collect  data  at  an  altitude,  Al.  Then  collect  data  at  a  lower  altitude, 
A2.  Though  the  aircraft  flies  at  different  altitudes,  it  tries  to  maintain  the  same  orbit. 
This  helps  image  registration  accuracy.  In  general,  two-pass  IFSAR  collection  and 
processing  is  not  as  easy  as  single-pass  IFSAR  collection  and  processing.  Two-pass 
IFSAR  collection  is  considered  to  be  useful  for  change  detection  (CD)  application  and 
measuring  terrain  motion  [9]. 
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Figure  2.12:  IFSAR  data  collection  geometry  with  single-pass.  In  this  collection, 

two  receivers  are  separated  vertically  on  the  same  platform  [9]. 


Figure  2.13:  IFSAR  data  collection  geometry  with  two-pass  interferometry.  In  this 
case,  the  radar  first  collects  data  at  an  altitude  Al.  Then  data  is  collected  at  a  lower 
altitude  A2.  In  both  altitudes,  the  platform  maintains  the  same  orbit  [9]. 
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2. 6  Summary 

This  section  surveys  SAR  signal  and  image  processing  concepts.  The  basic  radar 
operating  concept  was  presented.  Radar  signal  processing  using  a  radar  data  cube 
has  been  illustrated.  SAR  terminology  and  data  collection  geometry  have  been  dis¬ 
cussed.  SAR  range  and  cross-range  imaging  have  been  developed.  Finally,  a  pictorial 
description  of  SAR  data  collection  modalities  such  as  stripmap  SAR,  spotlight  SAR 
and  interferometric  SAR  (IFSAR)  is  illustrated. 
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III.  Circular  Synthetic  Aperture  Radar  (CSAR) 

This  chapter  describes  circular  synthetic  aperture  radar  (CSAR)  data  collection 
and  processing  concepts.  Benefits  of  CSAR  data  collection  have  been  discussed. 
A  very  important  CSAR  image  formation  algorithm  is  the  backprojection  algorithm. 
Principally,  two  different  backprojection  algorithms  can  be  found  in  the  literature; 
one  is  called  convolution  backprojection  (CBP)  [2,  5,  6,  20]  and  the  other  is  called 
the  backprojection  filtering  (BpjF)  [22].  In  this  chapter,  the  CBP  algorithm  will  be 
presented. 

3.1  CSAR  Processing 

SAR  data  collection  modalities  such  as  stripmap  SAR,  spotlight  SAR  and  in¬ 
terferometric  SAR  (IFSAR)  were  described  in  the  previous  chapter  and  use  a  radar 
trajectory  that  is  a  straight  flight  path.  These  data  collection  modalities  are  known 
as  linear  SAR  [19].  SAR  data  collection  using  a  circular  flight  path  is  known  as 
circular  SAR  (CSAR). 

In  CSAR,  the  radar  moves  around  a  circular  flight  path  while  steering  its  main- 
beam  to  constantly  illuminate  a  specific  ground  patch.  Antenna  steering  is  accom¬ 
plished  either  mechanically  or  electronically.  CSAR  data  collection  can  be  viewed  as 
spotlight  SAR  data  collection.  However,  the  difference  is  that  in  CSAR  the  data  is 
collected  using  a  circular  path  while  in  spotlight  SAR  data  is  collected  using  a  straight 
flight  path. 

CSAR  data  collection  geometry  is  presented  in  Figure  3.1.  This  circular  path 
span  is  from  0  to  360°  which  is  known  as  azimuth  coverage  or  aspect  angle,  (f).  The 
line  of  sight  (LOS)  distance  from  the  radar  to  the  imaging  scene  center  is  the  slant 
range.  The  angle  between  the  slant  range  to  the  ground  plane  is  called  the  elevation 
or  depression  angle,  6. 
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Figure  3.1:  Circular  SAR  (CSAR)  data  collection  geometry.  The  aircraft  carrying 
the  radar  flies  around  a  scene  of  interest.  The  radar  radiation  pattern  remains  constant 
on  the  scene. 

3.2  Merits  of  CSAR  Processing 

In  linear  SAR,  an  image  is  formed  using  data  collected  over  limited  azimuth  look 
angles.  This  imposes  two  problems  for  ATR  algorithm  development.  First,  images 
may  not  be  of  sufficient  quality  to  permit  target  identification  with  a  high  degree  of 
confidence  (e.g.  identification  of  a  school  bus  vs.  military  vehicle).  Second,  depending 
on  the  angular  orientation  of  the  target  in  the  scene,  it  will  be  difficult  to  determine 
the  target  type  (e.g.  identification  of  various  military  vehicle  of  similar  types)  [19]. 

CSAR  processing  provides  one  solution  to  the  aforementioned  problems.  First, 
it  is  well  understood  that  SAR  azimuth  resolution  (cross-range  resolution)  is  depen¬ 
dent  on  the  length  of  synthetic  aperture  [19].  Larger  apertures  improve  cross-range 
resolution.  CSAR  enables  large  aperture  sizes  by  steering  the  antenna  for  a  long 
period  and  therefore  provides  finer  azimuth  resolution  than  linear  SAR  [4], 

CSAR  has  the  advantage  of  imaging  the  scene  from  multiple  viewing  angles. 
This  flexibility  is  very  powerful  because  many  target  scattering  responses  depend  on 
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the  observation  angle.  Hence,  CSAR  imaging  from  multiple  viewing  angles  could 
provide  target  detection  and  identification  capability  with  high  degree  of  confidence. 

3.3  CSAR  Image  Formation  Techniques 

The  two  most  common  algorithms  for  CSAR  image  processing  are  the  polar 
format  algorithm  (PFA)  and  backprojection.  Data  collected  by  CSAR  systems  is 
called  Fourier  space  or  K-space  data  and  resides  on  an  annular  grid  as  shown  in 
Figure  3.2  [5]. 


Figure  3.2:  CSAR  data  representation  in  Fourier  domain  or  K-space  [5]. 

K-space  data  is  related  to  image  data  through  the  Fourier  transform.  Therefore, 
to  form  an  image,  2D  FFT  is  applied  to  the  K-space  data.  However,  before  apply¬ 
ing  the  2D  FFT,  an  FFT  based  imaging  method  (such  as  PFA)  needs  to  inscribe 
a  rectangular  grid  within  the  annular  region  as  shown  in  Figure  3.3.  This  is  called 
polar-to-rectangular  (Cartesian)  coordinate  conversion  of  K-space  data  and  is  accom¬ 
plished  by  applying  a  2D  interpolation  on  the  K-space  data.  Two  major  drawbacks  of 
the  PFA  includes:  (1)  all  radar  data  must  be  available  from  the  sensed  scene  to  apply 
polar  to  rectangular  conversion,  and  (2)  interpolation  inaccuracy  leads  to  degraded 
image  quality. 
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Figure  3.3:  Polar  to  rectangular  coordinate  interpolation  [5]. 

Backprojection  image  formation  uses  K-space  data  directly  and  does  not  require 
polar  to  rectangular  conversion.  Also,  all  radar  data  need  not  be  available  from  the 
sensed  scene  to  start  the  image  formation  process.  Further,  the  quality  of  images 
formed  by  backprojection  is  superior  to  that  of  the  PFA.  ffence,  images  formed  using 
backprojection  are  considered  the  Gold  Standard  for  CSAR  processing. 

3.4  Backprojection  Image  Formation 

The  backprojection  algorithm  was  conceived  in  the  field  of  computer-aided  to¬ 
mography  (CAT).  Later,  it  was  found  that  SAR  image  reconstruction  is  mathemat¬ 
ically  similar  to  CAT  and  the  researchers  developed  backprojection  for  SAR  [5]. 
Two  variations  of  backprojection  are  found  in  literature:  convolution  backprojection 
(CBP)  and  backprojection  filtering  (BpjF)  [22],  Some  papers  refer  CBP  as  filtered 
backprojection  (FBP)  or  factored  backprojection  (FBP)  [20]. 

In  CBP,  convolution  with  a  kernel  function  is  accomplished  on  each  projection. 
The  result  is  then  backprojected  to  a  grid.  In  the  BpjF  method,  first  backproject  the 
projections  on  a  image  grid.  Then  the  backprojected  image  is  filtered.  This  filtering 
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process  in  BpjF  introduces  image  reconstruction  inaccuracy  [22],  Hence  CBP  is  the 
preferred  method  for  backprojection. 

The  Radon  transform  is  fundamental  to  backprojection.  The  Radon  transform 
is  described  here  as  excerpted  from  the  book  by  Jain  [8].  Figure  3.4  illustrates  the 
Radon  transform  concept.  The  Radon  transform  maps  spatial  domain  (x,  y )  to  the 
(s,  9)  domain  [8].  Each  point  in  the  (s,  6)  space  is  a  line  in  the  spatial  domain  (x,  y ). 
The  relationships  between  Cartesian  (x,  y),  Polar (r,  </>)  and  Radon  (s,  6)  coordinates 
are  given  by  the  following  equations: 


X  = 

r  cos  <p 

(3.1) 

y  = 

r  sin  (p 

(3.2) 

s  = 

r  cos {9  —  <p) 

(3.3) 

Figure  3.4:  Radon  transform  concept.  Radon  transform  maps  spatial  domain  (x,y) 
or  polar  coordinate  (r,  (p )  into  the  (s,  9)  space  [8]. 


From  the  Radon  transform,  backprojection  is  defined  by  the  following  equations 
in  spatial  domain  and  polar  coordinate  respectively: 
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f{x,y) 


o 


(3.4) 


7T 


g(x  cos  9  +  yam  9,  6)d0 


f(r,<f>)  =  I  g(rcos(9-  <j>),6)d6 
o 


(3.5) 


Figure  3.5:  Backprojection  at  pixel  location  (x,  y )  or  (r,  <j) )  from  two  projections 
gi  and  g2  [8]. 

The  analytic  expressions  above  indicate  that  backprojection  is  the  summation 
or  integration  of  rays  that  pass  through  a  point  (x,  y)  or  (r,  </>).  In  SAR,  each  ray 
represents  a  range  profile  generated  from  a  radar  received  signal  or  pulse.  Figure  3.5 
shows  the  backprojection  concept  using  two  projections  g±  and  g2. 

Consider,  g(s,9)  =  gi(s)  +  g2(s).  Then  the  backprojection  at  (r,0)  is  given  by: 


=  Jl(Sl)  +  J2<S2) 


(3.6) 


where  si  =  rcos(0i  —  </>)  and  s2  =  rcos(d2  —  </>) 

Figure  3.6  is  a  visual  depiction  of  the  backprojection  image  formation  process. 
As  discussed  earlier,  rays  or  pulses  from  all  different  directions  are  backprojected 
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onto  a  rectangular  grid.  Summation  of  the  intensity  on  each  grid  is  the  resulting  pixel 
value. 

A 

o  O 

Object,  f{x,y) 

Figure  3.6:  Backprojection  image  reconstruction.  Projections  from  all  directions 

are  accumulated  to  form  an  image  of  the  objects  [8]. 

Having  presented  the  backprojection  concept  from  the  view  of  a  Radon  trans¬ 
form,  the  following  is  an  analytic  description  of  synthetic  aperture  radar  (SAR)  con¬ 
volution  backprojection  (CBP).  CBP  can  be  derived  from  the  2D  Fourier  transform. 
For  a  scene  with  ground  reflectivity  of  g(x,y),  the  2D  Fourier  transform  is  defined  by 
the  following  equations  [5]: 


G(u,  v ) 

g(x,y) 


(2ir) 


g(x,i i)e-^ux  +  '”l)dxdy 

G(u,v)esi'“+'v>dudv 


(3.7) 

(3.8) 


where  S  is  the  support  region  of  the  function  G(u,  v )  and  T  is  the  support  region  of 
the  function  g(x,  y ).  Consider  (a;  ,6)  is  the  representation  of  polar  coordinates  on  the 
(u,  v )  grid  shown  in  Figure  3.3.  Then  G(u  ,6  )  represents  the  value  of  G(u,  v )  along 
a  line  with  an  angle  9  along  the  u-axis. 
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Hence,  reconstructed  image  g^p  on  the  polar  coordinate  is  given  by  the  following 
equation  [5]: 

e/2  uj2 

gbp(x,y)=  I  J  G(u,  e)e~jw{-xcosd  +  y  sine)\cu\duj  dd  (3.9) 

-e/2  ui 


where  u  is  the  radar  frequency  and  6  is  the  projection  angle. 

From  the  above  equation,  SAR  image  reconstruction  using  CBP  can  be  viewed 
as  two  steps  process.  First,  each  radial  line  at  an  angle  6  is  multiplied  by  the  Jacobian, 
|w|.  This  is  the  convolution  step  since  multiplication  in  Fourier  space  is  equivalent  to 
convolution  in  time  domain.  This  result  is  then  backprojected  onto  an  image  grid. 


3. 5  Summary 

In  this  chapter  the  circular  synthetic  aperture  radar  (CSAR)  data  collection 
concept  has  been  presented.  Benefits  of  CSAR  processing  have  been  discussed.  Fi¬ 
nally,  an  important  CSAR  image  formation  technique,  the  backprojection  algorithm 
concept,  has  been  illustrated. 


IV.  Data  Generation  Methodology 


This  chapter  describes  the  data  generation  methodology  used  to  characterize  PSF 
responses  of  radially  displaced  scatterers.  The  SAR  data  collection  scenario 
simulated  for  this  research  is  the  CSAR  collection  method.  Data  is  collected  for 
full  azimuth  coverage  of  0  to  360°  by  placing  a  point  target  at  the  imaging  scene 
center  (0,0,0).  Similarly,  data  is  collected  by  placing  a  point  target  at  a  radially 
displaced  location  from  the  imaging  scene  center.  In  the  synthetic  data  generation 
process,  motion  compensation  which  is  a  phase  correction/adjustment  of  collected 
data  was  unnecessary.  However,  motion  compensation  is  essential  for  airborne  SAR 
data  collection  systems  given  that  aircraft  motion  distorts  signal  phase.  Details  of 
PSF  characteristics,  simulated  PSF  generation,  CSAR  data  collection  and  geometry, 
and  radar  waveform  selection  are  addressed  in  the  following  sections. 

4-1  Analytic  PSF  Characteristics 

Soumekh  [19]  provides  an  analytic  expression  for  the  2D  PSF  of  a  point  target 
located  at  the  origin,  i.e.  ( xn,yn )  =  (0,0).  The  2D  PSF  is  given  by: 

Table  4.1:  Experimental  parameters  used  to  characterize  PSF  of  a  radially  displaced 
point  scatterer.  In  this  table,  dr  represents  the  radial  displacements  of  the  point  target 
along  the  positive  x-axis  in  the  scene,  h  represents  location  of  the  image  plane  above 
the  ground  plane  where  h  —  0  is  the  ground  plane  and  9Z  is  the  depression  angle. 


Parameters 

Experimental  Range 

dr  (radial  displacement) 
h  (image  plane  height) 

9Z  (depression  angle) 

0,  3000,  5000  m 

0,  1,  2,  3  m 

15°,  30°,  45° 
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fn(x,y ) 

r 


Pn 


J\iX  Pmax)  Jliy  Pmin) 

Pmin 


Vx2  +  y2  , 


(4.1) 

(4.2) 


Pmax 


^kmax  COS  ($2)  5 


(4.3) 


Pmin 


^^min  COS  {0 


(4.4) 


where  J\  is  a  first-order  Bessel  function  of  the  first  kind,  r  is  the  scene  radius,  kmax 
is  maximum  wave  number  of  radar  signal,  kmin  is  minimum  wave  number  of  radar 
signal,  and  9Z  is  average  depression  angle  of  the  target  area. 

The  3D  PSF  of  a  target  located  at  the  scene  center  (0,0,0)  is  given  by  the 
following  equation  [13]: 

/(*,**)=^ (4.5) 

where  k\  and  k2  are  dehned  by  respectively  and  /i,  /2  denote  minimum  and 

maximum  radar  operating  frequencies. 

The  analytic  expressions  in  (4.1)  and  (4.5)  show  that  point  target  PSF  behavior 
is  a  function  of  radar  frequency  and  depression  angle.  Since  analytic  PSF  expressions 
for  radially  displaced  point  targets  are  difficult  to  obtain  [19],  this  study  focuses  on 
PSF  characterization  via  computer  simulation. 

The  PSF  of  a  radially  displaced  point  scatterer  can  be  characterized  as 

psf  =  f(dr,h,9z)  (4.6) 

where  dr  is  the  radial  displacement  of  the  point  target,  h  is  the  image  plane  height 
above  the  ground  and  9Z  is  the  depression  angle. 

Characterization  is  an  observation  of  changes  in  the  PSF  for  values  according 
to  the  Table  4.1  parameters.  For  the  on-center  target  (dr  =  0  m),  the  observation 
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Tabic  4.2:  Simulated  parameters  used  in  backprojection  imaging.  A  bandwidth  Aw 
of  0.5  GHz  provides  a  range  resolution  Ar  of  0.3  m,  a  center  frequency  of  /o  2.5  GHz 
provides  a  wavelength  A  of  0.12  m. 


Parameters 

Value 

Acu  (bandwidth) 

0.5  GHz 

/o  (center  frequency) 

2.5  GHz 

A  (wavelength) 

0.12  nr 

r  (scene  radius) 

7000  nr 

A  (sensor  altitude) 

7000  nr 

9Z  (depression  angle) 

45° 

Ar.  (range  resolution) 

0.3  nr 

includes  the  PSF  response  is  circular  and  symmetric  about  the  2-axis  or  image  plane 
height  h.  For  the  radially  displaced  target  (dr  =  3000  m  and  5000  m),  the  observation 
includes  the  PSF  response  is  elliptic  and  asymmetric  about  the  2-axis  or  image  plane 
height  h. 

4-2  Simulated  PSF  Generation 

Table  4.2  provides  the  simulation  parameters.  Parameters  are  chosen  for  conve¬ 
nience.  However,  other  parameters  can  be  chosen  as  well.  Simulations  were  conducted 
by  placing  a  point  scatterer  at  the  imaging  scene  center  and  forming  an  image  (PSF) 
of  the  scatterer  using  a  backprojection  imaging  algorithm.  The  simulations  were  then 
run  using  a  radially  displaced  point  scatterer,  i.e. ,  a  point  scatterer  displaced  from 
the  imaging  center  along  the  positive  cc-axis  by  a  predetermined  distance.  Simulations 
were  run  for  these  two  cases  using  a  depression  angle  of  6Z  =  45°,  bandwidth  of  Ac o  = 
0.5  GHz  and  azimuth  coverage  of  0  to  360°.  A  bandwidth  of  Ac o  =  0.5  GHz  provides 
a  range  resolution  of  Ar  «  0.3  meters.  A  center  frequency  of  /o  =  2.5  GHz  was  used 
and  corresponds  to  a  wavelength  of  A  =  0.12  m. 

Backprojection  image  formation  was  accomplished  by  coherently  summing  all 
collected  returns  for  a  given  depression  angle  using  a  full  azimuth  span.  Images 
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formed  using  the  full  360°  azimuth  span  provide  the  greatest  amount  of  target  infor¬ 
mation/image  detail  because  the  target  has  been  observed  from  all  aspect  angles. 

4-3  Data  Collection  Geometry 

z 


Figure  4.1:  CSAR  collection  and  imaging  geometry.  CSAR  equipped  aircraft  flies 
in  a  circular  path  of  radius  7000  m  at  an  altitude  of  7000  m  from  the  ground  plane. 

The  simulated  CSAR  imaging  geometry  is  shown  in  Fig.  4.1.  As  indicated,  the 
CSAR  equipped  aircraft  flies  a  circular  path  around  the  2-axis  at  a  depression  angle 
9Z  measured  from  the  X-Y  plane.  For  this  data  collection,  depression  angle  6Z  is  set 
to  45°.  The  aircraft  flies  in  a  circular  path  of  radius  r  =  7000  m  at  an  altitude  A  = 
7000  m  from  the  ground  plane.  Slant  range  R  from  the  image  scene  center  (0,0,0)  to 
the  aircraft  is  approximately  9900  m. 

Data  was  collected  over  an  azimuth  span  of  0  to  360°  under  two  conditions, 
including:  (1)  a  point  target  at  the  origin  (0,0,0)  and  (2)  a  radially  displaced  point 
target  located  at  3000  m  and  5000  m  along  the  positive  x-axis. 


32 


Figure  4.2:  Sine  shaped  compressed  radar  waveform. 

4-4  Radar  Waveform 

Time  domain  signal  processing  was  used  for  all  data  generated  in  support  of  this 
study.  In  other  words,  frequency  domain  signal  representation  and  comprehensive 
waveform  design,  analysis,  and  pulse  compression  methods  are  not  employed.  The 
radar  antenna  parameters,  beamwidth  and  polarization,  have  not  been  considered  for 
the  simulation. 


Figure  4.3:  Gaussian  shaped  compressed  radar  waveform. 
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Figure  4.4:  CSAR  image  formed  using  a  sine  shaped  compressed  radar  waveform. 

For  an  azimuth  span  of  0  to  360° ,  one  radar  pulse  was  transmitted  every  0.5°. 
Thus  720  pulses  were  transmitted  and  received  by  the  radar.  Two  different  compressed 
waveforms  have  been  used  to  study  radar  waveform  impact  on  target  signature:  one 
is  a  sine  shaped  waveform  and  the  other  is  the  Gaussian  shaped  waveform.  Figure  4.2 
and  Fig.  4.3  show  sine  and  Gaussian  shaped  waveforms,  respectively.  Figure  4.4  and 
Fig.  4.5  show  two  CSAR  images  formed  using  each  of  the  aforementioned  waveforms. 
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Figure  4.5:  CSAR  image  formed  using  a  gaussian  shaped  compressed  radar  wave¬ 
form. 
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Since  a  sine  shaped  waveform  represents  a  more  realistic  radar  waveform  model  for  a 
practical  radar  system,  the  Gaussian  shaped  waveform  is  not  used  for  data  generation. 

4-5  Summary 

This  chapter  outlines  the  data  generation  method  for  PSF  characterization. 
Data  is  collected  over  an  azimuth  span  of  0  to  360°  by  placing  a  point  scatterer  at 
the  imaging  scene  center  (0,0,0)  as  well  as,  by  placing  the  point  scatterer  in  a  radially 
displaced  along  the  positive  x-axis.  Key  CSAR  parameters  for  this  simulation  include 
a  bandwidth  of  Ac o  =  0.5  GHz,  center  radar  frequency  of  /o  =  2.5  GHz,  and  a 
depression  angle  of  6Z  =  45°.  The  aircraft  flies  a  circular  path  of  radius  7000  m  at  an 
altitude  of  7000  nr  from  the  ground  plane.  The  radar  uses  a  sine  shaped  time  domain 
compressed  waveform. 
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V.  Results  and  Analysis 


This  chapter  presents  results  and  analysis  describing  the  point  spread  function 
(PSF)  characteristics  of  a  point  scatterer.  First,  key  observations  characterizing 
the  PSF  of  a  radially  displaced  scatterer  is  presented.  Then,  baseline  PSF  results  for 
a  point  target  located  at  the  scene  center  (0,0,0)  are  demonstrated.  Both  2D  PSFs 
and  3D  plots  of  the  PSFs  for  the  on-center  point  target  are  presented  and  analyzed. 
PSF  results  for  radially  displaced,  off-center  point  targets,  are  then  presented.  Again, 
both  2D  PSFs  and  3D  plots  of  the  PSF  responses  are  illustrated  and  analyzed. 

5.1  PSF  Characteristics  of  an  Off-center  Target 

Table  5.1  depicts  key  observations  characterizing  the  PSF  of  an  off-center  point 
target.  With  a  radial  displacement  dr  of  3000  m  from  the  scene  center  (dr  =  0  m),  the 
PSF  of  the  point  target  tends  to  become  an  elliptic  shape  from  a  circular  shape.  With 
an  additional  displacement  of  dr  =  5000  m,  the  ellipticity  increases.  Therefore,  as  A dr 
increases,  the  PSF  becomes  more  elliptic  shaped  and  asymmetric  about  the  z-axis.  A 
change  in  the  image  plane  from  h  =  0  m  to  h  =  3  m,  the  PSF  response  changes  from  a 
fully  focused  image  to  a  degraded  annulus  shaped  image.  When  h  —  0  in  (i.e.  ground 
plane  is  the  image  plane),  the  PSF  is  fully  focused.  As  Ah  increases,  the  width  of 
the  annulus  decreases.  This  is  an  indication  that  a  target  is  present  but  it’s  image 
is  not  focused.  Also  for  h  >  0,  the  peak  PSF  response  is  no  longer  observed  at  the 
true  target  location.  A  change  in  9Z,  changes  the  radius  of  the  annulus.  Maximum 
increase  in  the  annulus  radius  occurs  at  9Z  =  45°. 

5.2  PSF  Results:  Point  target  at  the  origin 

The  PSF  of  an  on-center  (0,0,0)  point  target  serves  as  a  basis  for  characterizing 
PSF  of  displaced  point  targets.  Hence,  this  section  provides  a  detailed  analysis  of  the 
on-center  point  target  PSF. 

First,  the  2D  PSF  of  a  scatterer  in  the  ground  plane  is  presented.  Then,  2D 
PSFs  in  different  image  planes  at  heights  1,  2,  and  3  meters  from  the  ground  plane 
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Table  5.1:  Results  characterizing  PSF  of  a  radially  displaced  point  scatterer.  In  this 
table,  A dr  represents  changes  in  the  radial  displacements  of  the  point  target  along 
the  positive  x-axis  in  the  scene,  Ah  represents  changes  in  the  locations  of  the  image 
plane  above  the  ground  plane  where  h  —  0  is  the  ground  plane  and  A 9Z  represents 
changes  in  the  depression  angle. 


Parameters 

Observation 

A dr  (radial  displacement) 

PSF  changes  from  a  circle  to  an  ellipse  and  PSF 
becomes  asymmetric  about  the  z-axis  as  dr  increases. 

Ah  (image  plane) 

Width  of  the  annulus  decreases  as  h  increases. 

Also  for  h  >  0,  the  peak  PSF  response  is  no  longer 
observed  at  the  true  target  location. 

A 6Z  (depression  angle) 

Annulus  radius  changes.  Maximum  occurs  at  9Z  =  45° 

are  presented.  Similarly,  the  3D  plot  of  the  PSF  of  a  scatterer  at  the  ground  plane  is 
presented.  Then,  3D  plots  of  the  PSFs  in  different  image  planes  at  heights  1,  2,  and 
3  meters  from  the  ground  plane  are  presented. 


x  (meter) 


Figure  5.1:  2D  PSF  for  point  scatterer  located  at  the  scene  center  (0,0,0).  The 

image  was  formed  for  a  depression  angle  6Z  of  45°,  a  bandwidth  Ac j  of  0.5  GHz  and 
azimuth  span  of  0  to  360°  [12]. 
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5.2.1  2D  PSF  Responses.  Figure  5.1  through  Fig.  5.5  show  baseline  PSF 
results  for  a  point  target  located  at  the  origin  (0,0,0)  or  radial  displacement  dr  =  0  m. 
Figure  5.1  shows  the  2D  PSF  of  a  point  scatterer  located  at  the  scene  center  using 
bandwidth  Ac o  of  0.5  GHz,  and  full  azimuth  coverage  of  0  to  360°.  The  image  plane 
is  the  ground  plane. 

Figure  5.2  through  Fig.  5.4  show  2D  PSFs  formed  at  different  heights  above  the 
ground  plane  for  a  point  scatterer  located  at  the  scene  center  (0,0,0).  The  images 
were  formed  for  a  depression  angle  of  45°  and  full  azimuth  coverage  of  0  to  360°.  PSF 
results  at  height  h  —  1,  2  and  3  meters  are  presented  in  Fig.  5.2  through  Fig.  5.4 
respectively.  As  shown,  for  all  focusing  heights  ( h )  greater  than  zero  (i.e.  above 
the  ground  plane),  the  response  degrades  to  an  annulus  shape.  The  maximum  PSF 
response  within  the  annulus  decreases  with  h  and  are  all  20  dB  or  more  below  the 
peak  h  =  0  ground  plane  response  in  Fig.  5.1.  The  annulus  radius  is  a  function  of 
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Figure  5.2:  2D  PSF  for  the  point  scatterer  located  at  the  scene  center  (0,0,0).  The 
image  plane  is  1  meter  above  the  ground  plane.  The  image  was  formed  for  a  depression 
angle  6L  of  45°,  a  bandwidth  Aou  of  0.5  GHz  and  azimuth  span  of  0  to  360°  [12], 
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depression  angle  (Qz)  and  height  (h)  above  the  point  target  and  is  given  by  [12]: 


r  =  h  ■  tan(#z) 


(5.1) 


Figure  5.5  shows  the  X-Z  cut  of  the  3D  PSF  results  with  the  point  scatterer 
located  at  the  scene  center.  This  figure  illustrates  the  fact  that  for  an  on-center 
target  (i.e.  when  target’s  radial  displacement  dr  =  0),  the  PSF  is  symmetric  about 
the  z-axis. 
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Figure  5.3:  2D  PSF  for  the  point  scatterer  located  at  the  scene  center  (0,0,0). 

The  image  plane  is  2  meters  above  the  ground  plane.  The  image  was  formed  for  a 
depression  angle  6~  of  45°,  a  bandwidth  Acn  of  0.5  GHz  and  azimuth  span  of  0  to 
360°  [12]. 
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Figure  5.4:  2D  PSF  for  the  point  scatterer  located  at  the  scene  center  (0,0,0). 

The  image  plane  is  3  meters  above  the  ground  plane.  The  image  was  formed  for  a 
depression  angle  9Z  of  45°,  a  bandwidth  Au  of  0.5  GHz  and  azimuth  span  of  0  to 
360°  [12]. 
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Figure  5.5:  X-Z  plane  cut  of  the  3D  PSF  for  point  scatterer  located  at  the  scene 
center  (0,0,0).  The  image  was  formed  for  a  depression  angle  9Z  of  45°,  a  bandwidth 
Au  of  0.5  GHz  and  azimuth  span  of  0  to  360°  [12]. 
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5.2.2  3D  Plots  of  the  PSF  Responses.  Figure  5.6  through  Fig.  5.9  show 
3D  plots  of  the  PSFs  formed  at  different  heights  above  the  ground  plane  for  a  point 
scatterer  located  at  the  scene  center  (0,0,0)  or  radial  displacement  dr  =  0  m.  The 
images  were  formed  for  a  depression  angle  6Z  of  45°  and  full  azimuth  coverage  of  0  to 
360°.  PSF  results  at  height  h  =  0,  1,  2  and  3  meters  are  presented  in  Fig.  5.6  through 
Fig.  5.9  respectively,  where  h  —  0  is  the  ground  plane.  As  shown,  when  focusing 
at  heights  greater  than  zero,  the  response  degrades  to  an  annulus  shape.  Fig.  5.10 
shows  the  3D  PSF  results  for  CSAR  conditions  with  the  point  scatterer  located  at 
the  scene  center.  Consistent  with  published  results  [1,13],  the  3D  PSF  is  cone  shaped 
and  symmetric  about  the  z-axis.  For  all  h,  the  peak  PSF  response  remains  at  the 
true  target  location. 
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Figure  5.6:  3D  plot  of  PSF  for  the  point  scatterer  located  at  the  scene  center  (0,0,0). 
The  image  plane  is  the  ground  plane.  The  image  was  formed  for  a  depression  angle 
9Z  of  45°,  a  bandwidth  A u  of  0.5  GHz  and  azimuth  span  of  0  to  360° 
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Figure  5.7:  3D  plot  of  PSF  for  the  point  scatterer  located  at  the  scene  center 

(0,0,0).  The  image  plane  is  1  meter  above  the  ground  plane.  The  image  was  formed 
for  a  depression  angle  6Z  of  45°,  a  bandwidth  Au>  of  0.5  GHz  and  azimuth  span  of  0 
to  360° 
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Figure  5.8:  3D  plot  of  PSF  for  the  point  scatterer  located  at  the  scene  center  (0,0,0). 
The  image  plane  is  2  meters  above  the  ground  plane.  The  image  was  formed  for  a 
depression  angle  9Z  of  45°,  a  bandwidth  Ac v  of  0.5  GHz  and  azimuth  span  of  0  to  360° 
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Figure  5.9:  3D  plot  of  PSF  for  the  point  scatterer  located  at  the  scene  center  (0,0,0). 
The  image  plane  is  3  meters  above  the  ground  plane.  The  image  was  formed  for  a 
depression  angle  9Z  of  45°,  a  bandwidth  Acn  of  0.5  GHz  and  azimuth  span  of  0  to  360° 
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Figure  5.10:  3D  PSF  for  point  scatterer  located  at  the  scene  center  (0,0,0).  The 
image  was  formed  for  a  depression  angle  9Z  of  45°,  a  bandwidth  Acj  of  0.5  GHz  and 
azimuth  span  of  0  to  360° 
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5.3  PSF  Results:  Radially  displaced  point  target 

As  mentioned  earlier  in  Chapter  4,  an  analytic  closed-form  solution  describing 
the  PSF  of  a  radially  displaced  target  is  difficult  to  obtain  [19].  Hence,  PSF  results 
presented  here  for  a  radially  displaced  point  target  is  based  on  computer  simulation. 
First,  the  2D  PSF  of  a  scatterer  at  the  ground  plane  (i.e.  h  —  0  m)  is  presented.  Then, 
2D  PSFs  at  different  image  planes  located  at  heights  of  1,  2,  and  3  meters  above  the 
ground  plane  are  presented.  Similarly,  the  3D  plot  of  the  PSF  of  a  scatterer  at  the 
ground  plane  is  presented.  Then,  3D  plots  of  PSFs  at  different  image  planes  located 
at  heights  1,  2,  and  3  meters  above  the  ground  plane  is  presented. 

5.3.1  2D  PSF  Responses.  Figure  5.11  through  Fig.  5.14  show  2D  PSF  results 
for  point  target  displaced  dr  =  3000  meters  along  the  positive  x-axis.  Figure  5.11 
shows  the  2D  PSF  of  a  point  scatterer  using  a  depression  angle  6Z  of  45°,  a  bandwidth 
Au;  of  0.5  GHz,  and  full  azimuth  coverage  of  0  to  360°.  The  image  plane  is  the  ground 
plane. 

Fig.  5.12  through  Fig.  5.14  show  2D  PSFs  formed  at  different  heights  above 
the  ground  plane  for  a  point  scatterer  displaced  3000  meters  along  the  positive  x- 
axis.  The  images  were  formed  for  a  depression  angle  of  9Z  =  45°  and  full  azimuth 
coverage  of  0  to  360°.  In  these  cases,  the  previous  annulus  responses  become  elliptical 
and  asymmetric  about  the  z-axis.  PSF  results  at  height  h  =  1,  2,  and  3  meters  are 
presented  in  Fig.  12  through  Fig.  14,  respectively.  For  h  >  0,  the  peak  response  is  no 
longer  observed  at  the  true  target  location  (i.e.  dr  =  3000  m).  Figure  5.15  shows  the 
X-Z  cut  of  the  3D  PSF  results  for  the  point  scatterer  displaced  3000  meters  along 
the  positive  x-axis.  In  this  case,  previous  cone  shown  in  Fig.  5.10  becomes  a  tilted 
cone  and  thus  X-Z  cut  is  tilted. 
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Figure  5.11:  2D  PSF  for  point  scatterer  displaced  3000  meters  along  the  positive 
rc-axis.  The  image  plane  is  the  ground  plane.  The  image  was  formed  for  a  depression 
angle  9~  of  45°,  a  bandwidth  Ac j  of  0.5  GHz  and  azimuth  span  of  0  to  360°  [12]. 
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Figure  5.12:  2D  PSF  for  point  scatterer  displaced  3000  meters  along  the  positive 
cr-axis.  The  image  plane  is  1  meter  above  the  ground  plane.  The  image  was  formed 
for  a  depression  angle  6Z  of  45°,  a  bandwidth  Acu  of  0.5  GHz  and  azimuth  span  of  0 
to  360°  [12]. 
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Figure  5.14:  2D  PSF  for  point  scatterer  displaced  3000  meters  along  the  positive 
x-axis.  The  image  plane  is  3  meters  above  the  ground  plane.  The  image  was  formed 
for  a  depression  angle  6h  of  45°,  a  bandwidth  Acu  of  0.5  GHz  and  azimuth  span  of  0 
to  360°  [12]. 


Figure  5.13:  2D  PSF  for  point  scatterer  displaced  3000  meters  along  the  positive 
x-axis.  The  image  plane  is  2  meters  above  the  ground  plane.  The  image  was  formed 
for  a  depression  angle  6Z  of  45°,  a  bandwidth  Au>  of  0.5  GHz  and  azimuth  span  of  0 
to  360°  [12]. 
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Figure  5.15:  X-Z  plane  cut  of  the  3D  PSF  for  point  scatterer  displaced  3000  meters 
along  the  positive  x-axis.  The  image  was  formed  for  a  depression  angle  9Z  of  45°,  a 
bandwidth  A u  of  0.5  GHz  and  azimuth  span  of  0  to  360°  [12]. 


Figure  5.16  through  Fig.  5.19  show  2D  PSF  results  for  point  target  displaced 
dr  =  5000  m  along  the  positive  x-axis.  The  images  were  formed  for  a  depression 
angle  9Z  of  45°,  full  azimuth  coverage  of  0  to  360°  and  at  height  h  =  0,1,  2  and  3 
meters  respectively  from  the  ground  plane.  With  the  additional  radial  displacement, 
the  PSF  tends  to  become  more  elliptically  shaped.  For  h  >  0,  PSF  responses  degrade 
to  annulus  shapes  and  asymmetric  about  the  z-axis.  Also  for  h  >  0,  the  peak  PSF 
response  is  no  longer  observed  at  the  true  target  location  (i.e.  dr  =  5000  m). 
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Figure  5.16:  2D  PSF  for  point  scatterer  displaced  5000  meters  along  the  positive 
tr-axis.  The  image  plane  is  the  ground  plane.  The  image  was  formed  for  a  depression 
angle  9Z  of  45°,  a  bandwidth  Ac o  of  0.5  GHz  and  azimuth  span  of  0  to  360° 
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Figure  5.17:  2D  PSF  for  point  scatterer  displaced  5000  meters  along  the  positive 
tr-axis.  The  image  plane  is  1  meter  above  the  ground  plane.  The  image  was  formed 
for  a  depression  angle  9Z  of  45°,  a  bandwidth  Au>  of  0.5  GHz  and  azimuth  span  of  0 
to  360° 
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Figure  5.18:  2D  PSF  for  point  scatterer  displaced  5000  meters  along  the  positive 
x-axis.  The  image  plane  is  2  meters  above  the  ground  plane.  The  image  was  formed 
for  a  depression  angle  9Z  of  45°,  a  bandwidth  Au>  of  0.5  GHz  and  azimuth  span  of  0 
to  360° 
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Figure  5.19:  2D  PSF  for  point  scatterer  displaced  5000  meters  along  the  positive 
x-axis.  The  image  plane  is  3  meters  above  the  ground  plane.  The  image  was  formed 
for  a  depression  angle  9Z  of  45°,  a  bandwidth  Au>  of  0.5  GHz  and  azimuth  span  of  0 
to  360° 
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5.3.2  3D  plots  of  the  PSF  Responses.  Figure  5.20  through  Fig.  5.23  show  3D 
plots  of  PSFs  formed  at  different  heights  above  the  ground  plane  for  a  point  scatterer 
displaced  dr  =  3000  meters  along  the  positive  x-axis.  The  images  were  formed  for  a 
depression  angle  9Z  of  45°  and  full  azimuth  coverage  of  0  to  360°.  PSF  results  at  height 
h  =  0,  1,  2  and  3  meters  are  presented  in  Fig.  5.20  through  Fig.  5.23,  respectively, 
where  h  —  0  is  the  ground  plane.  As  shown,  when  focusing  heights  are  greater  than 
zero  (i.e.  image  plane  is  above  the  ground  plane),  the  previous  annulus  responses 
become  elliptically  shaped  and  asymmetric  about  the  z-axis.  Also  for  h  >  0,  3D  plots 
of  the  PSF  responses  become  tilted. 
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Figure  5.20:  3D  plot  of  PSF  for  the  point  scatterer  displaced  3000  meters  along 
the  positive  x-axis.  The  image  plane  is  the  ground  plane.  The  image  was  formed  for 
a  depression  angle  9Z  of  45°,  a  bandwidth  Ac v  of  0.5  GHz  and  azimuth  span  of  0  to 
360° 
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Figure  5.21:  3D  plot  of  PSF  for  the  point  scatterer  displaced  3000  meters  along 
the  positive  x-axis.  The  image  plane  is  1  meter  above  the  ground  plane.  The  image 
was  formed  for  a  depression  angle  9Z  of  45°,  a  bandwidth  A u>  of  0.5  GHz  and  azimuth 
span  of  0  to  360° 
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Figure  5.22:  3D  plot  of  PSF  for  the  point  scatterer  displaced  3000  meters  along 
the  positive  x-axis.  The  image  plane  is  2  meters  above  the  ground  plane.  The  image 
was  formed  for  a  depression  angle  9Z  of  45°,  a  bandwidth  A u>  of  0.5  GHz  and  azimuth 
span  of  0  to  360° 
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Figure  5.23:  3D  plot  of  PSF  for  the  point  scatterer  displaced  3000  meters  along 
the  positive  x-axis.  The  image  plane  is  3  meters  above  the  ground  plane.  The  image 
was  formed  for  a  depression  angle  9Z  of  45°,  a  bandwidth  Au>  of  0.5  GHz  and  azimuth 
span  of  0  to  360° 


Figure  5.24  through  Fig.  5.27  show  3D  plots  of  PSFs  formed  at  different  heights 
above  the  ground  plane  for  a  point  scatterer  displaced  dr  =  5000  meters  along  the 
positive  x-axis.  The  images  were  formed  for  a  depression  angle  of  9Z  =  45°  and  full 
azimuth  coverage  of  0  to  360°.  3D  PSF  results  at  height  h  —  0,  1,  2  and  3  meters 
are  presented  in  Fig.  5.24  through  Fig.  5.27  respectively,  where  h  =  0  is  the  ground 
plane.  With  the  additional  radial  displacement,  the  PSF  ellipticity  increases  relative 
to  previous  cases  and  asymmetric  about  the  x-axis.  3D  plot  of  the  PSF  response 
becomes  more  tilted. 
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Figure  5.24:  3D  plot  of  PSF  for  the  point  scatterer  displaced  5000  meters  along 
the  positive  x-axis.  The  image  plane  is  the  ground  plane.  The  image  was  formed  for 
a  depression  angle  9Z  of  45°,  a  bandwidth  Au  of  0.5  GHz  and  azimuth  span  of  0  to 
360° 
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Figure  5.25:  3D  plot  of  PSF  for  the  point  scatterer  displaced  5000  meters  along 
the  positive  x-axis.  The  image  plane  is  1  meter  above  the  ground  plane.  The  image 
was  formed  for  a  depression  angle  9Z  of  45°,  a  bandwidth  A  to  of  0.5  GHz  and  azimuth 
span  of  0  to  360° 
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Figure  5.26:  3D  plot  of  PSF  for  the  point  scatterer  displaced  5000  meters  along 
the  positive  x-axis.  The  image  plane  is  2  meters  above  the  ground  plane.  The  image 
was  formed  for  a  depression  angle  9Z  of  45°,  a  bandwidth  Au  of  0.5  GHz  and  azimuth 
span  of  0  to  360° 


Figure  5.27:  3D  plot  of  the  PSF  for  the  point  scatterer  displaced  5000  meters  along 
the  positive  x-axis.  The  image  plane  is  3  meters  above  the  ground  plane.  The  image 
was  formed  for  a  depression  angle  9~  of  45°,  a  bandwidth  A  to  of  0.5  GHz  and  azimuth 
span  of  0  to  360° 
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5.4  Summary 

This  chapter  provides  results  describing  the  PSF  characteristics  of  a  point  scat¬ 
tered  The  PSFs  were  formed  for  a  depression  angle  9Z  of  45°,  a  bandwidth  A ui  of 
0.5  GHz  and  azimuth  span  of  0  to  360°.  Both  2D  PSFs  and  3D  plots  of  the  PSFs  are 
presented  for  a  point  target  located  at  the  scene  center  (0,0,0)  or  radial  displacement 
dr  =  0  m.  Then,  2D  PSFs  and  3D  plot  of  the  PSFs  for  the  radially  displaced  (dr  = 
3000,  5000  m  )  point  target  have  been  presented.  Results  show  that  for  an  on-center 
(dr  =  0  m)  point  target,  the  2D  PSF  is  focused  and  circular.  The  3D  PSF  response  is 
cone  shaped  and  symmetric  about  the  z-axis.  This  result  is  consistent  with  previously 
published  research  [1,13]. 

For  a  radially  displaced  point  scatterer,  the  PSF  becomes  elliptic.  As  radial 
displacement  dr  increases  from  3000  m  to  5000  m,  the  PSF  becomes  more  elliptic 
shaped  and  asymmetric  about  the  z-axis.  A  change  in  the  image  plane  from  h  =  0  m 
to  h  =  3  m,  the  PSF  response  changes  from  a  fully  focused  image  to  a  degraded 
annulus  shaped  image.  Also  for  h  >  0,  the  peak  PSF  response  is  no  longer  observed 
at  the  true  target  location.  A  change  in  the  depression  angle  9Z,  changes  the  annulus 
radius.  The  maximum  increase  in  the  annulus  radius  occurs  when  0Z  =  45°. 
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VI.  Conclusions 


The  Point  Spread  Function  (PSF)  of  a  radially  displaced  point  scatterer  is  char¬ 
acterized  for  CSAR  imagery  via  simulation.  Backprojection  image  formation 
was  used  to  simulate  the  PSF  of  the  on-center  and  off-center  point  targets.  Impor¬ 
tant  radar  parameters  used  for  this  research  include  a  radar  center  frequency  /0  of 
2.5  GHz,  range  resolution  Ar  of  0.3  m  and  depression  angle  6Z  of  45°.  Contribution 
of  this  research  to  SAR  application  has  been  introduced.  Finally,  future  research  on 
this  topic  has  been  presented. 

6.1  Research  Summary 

Simulated  2D  and  3D  PSF  results  are  presented  for  a  point  scatterer  using 
Circular  Synthetic  Aperture  Radar  (CSAR)  processing.  When  the  scatterer  is  located 
at  the  imaging  scene  center,  with  full  360°  azimuth  coverage  and  a  depression  angle 
6Z  of  45°  the  scattering  response  is  circular  and  symmetrical  around  the  z-axis.  The 
3D  PSF  response  is  cone  shaped  [1,11,13]. 

When  the  point  target  is  radially  displaced,  the  PSF  response  is  elliptic  and  no 
longer  symmetric  about  the  z-axis.  In  this  case,  the  previous  3D  cone  becomes  more 
tilted  as  the  distance  from  the  scene  center  increases.  Thus,  a  2D  elevation  slice  of 
the  tilted  cone  gives  a  conic  section.  For  the  cases  simulated,  the  conic  section  was 
an  ellipse  having  increased  ellipticity  with  larger  radial  displacement. 

Another  PSF  characteristic  was  observed  by  forming  2D  PSF  responses  at 
heights  of  h  =  1,  2  and  3  meters  above  the  ground  plane.  It  was  observed  that 
the  maximum  PSF  response  within  the  annulus  decreases  with  h.  For  example,  the 
peak  response  at  h  —  1  was  20  dB  or  more  below  the  peak  ground  plane  response 
at  h  =  0.  As  a  result,  the  annulus  responses  appeared  wider  at  the  lower  heights. 
For  the  point  scatterers  displaced  from  the  scene  center,  a  2D  cut  at  a  height  above 
the  ground  plane  illustrates  the  combined  effect  of  the  elliptical  shape  and  annulus 
broadening  as  depicted  in  Figure  5.17.  Also  for  the  radially  displaced  scatterer  when 
h  >  0,  the  peak  PSF  response  is  no  longer  observed  at  the  true  target  location.  It  was 
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also  observed  that  a  change  in  the  depression  angle  9Z,  changes  the  annulus  radius. 
The  maximum  increase  in  the  annulus  radius  occurs  at  9Z  of  45°, 

6.2  Contribution  to  SAR  research 

As  discussed  in  the  introduction,  there  is  an  important  need  to  image  large  scene 
sizes  with  CSAR.  Unfortunately,  as  a  result  of  the  imaging  geometry,  the  PSF  is  not 
invariant  to  scatterer  position  in  the  scene.  In  order  to  perform  image  interpretation 
or  exploitation  of  large  scenes  using  CSAR,  it  is  important  to  understand  how  the 
PSF  changes  as  a  function  of  the  target’s  position  in  the  scene.  This  research  provides 
one  experimental  step  toward  obtaining  this  understanding. 

In  addition,  this  research  provides  insight  into  how  one  might  formulate  an 
analytic  solution  to  the  problem.  Further,  since  many  signal  and  image  processing 
approaches  to  detection  and  ATR  rely  on  knowledge  of  the  system  PSF,  results  of 
this  thesis  could  lead  to  new  and  improved  algorithms  for  CSAR  imaging  of  large 
scenes.  In  particular,  the  development  of  target  recognition  and  detection  approaches 
for  both  stationary  (e.g.  SAR)  and  moving  targets  (e.g.,  GMTI)  could  benefit  from 
this  research. 

6.3  Future  Research  on  Radially  Displaced  Targets  in  CSAR 

As  stated  earlier,  this  research  endeavor  characterized  the  PSF  of  a  radially 
displaced  target  in  CSAR  imagery.  The  following  section  highlights  potential  follow- 
on  research  to  build  on  the  fundamentals  of  this  thesis  topic. 

6.3.1  CSAR  response  to  canonical  shapes.  One  direction  of  research  could 
be  to  investigate  the  image  response  of  some  canonical  shapes  at  radially  displaced 
locations.  Canonical  shapes  such  as  dihedrals,  trihedrals,  cylinders  would  be  inter¬ 
esting  and  important  shapes  to  consider.  Understanding  these  responses  would  be 
an  important  step  towards  understanding  the  images  of  complex  targets  such  as  the 
buildings,  motor  vehicles,  etc.  leading  to  many  persistent  staring  radar  applications.. 
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6.3.2  Analytic  PSF  expression  for  radially  displaced  point  scatterers.  Re¬ 
search  has  been  conducted  to  find  an  analytic  closed-form  solution  describing  the  PSF 
of  a  point  target  centered  at  origin  [13].  This  solution  involves  a  formulation  that  in¬ 
cludes  Bessel  function  evaluation.  For  a  radially  displaced  point  target,  an  analytical 
closed-form  solution  has  not  yet  been  developed.  One  approach  towards  developing 
an  analytical  solution  could  involve  transforming  the  on-center  [13]  analytical  solu¬ 
tion  based  on  insight  provided  in  this  thesis  research.  Another  avenue  of  investigation 
could  include  the  development  of  a  numerical  approximation  that  describes  the  PSF 
of  a  radially  displaced  target. 

6.3.3  Detection  and  ATR  algorithms  based  on  PSF  response.  Many  signal 
processing  algorithms  take  advantage  of  a  priori  PSF  knowledge.  Algorithms  using 
this  PSF  knowledge  may  be  a  fruitful  area  of  research. 
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